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ABSTRACT 

The present paper describes the analysis of multiple RXTE/PCA data of the 
black hole binary with superluminal jet, XTE J1550 — 564, acquired during its 
1999-2000 outburst. The X-ray spectra show features typical of the high/soft 
spectral state, and can approximately be described by an optically thick disk 
spectrum plus a power-law tail. Three distinct spectral regimes, named standard 
regime, anomalous regime, and apparently standard regime, have been found 
from the entire set of the observed spectra. When the X-ray luminosity is well 
below ~ 6 x 10 38 erg s _1 (assuming a distance of 5 kpc), XTE J1550 — 564 
resides in the standard regime, where the soft spectral component dominates 
the power-law component and the observed disk inner radius is kept constant. 
When the luminosity exceeds the critical luminosity, the apparently standard 
regime is realized, where luminosity of the optically thick disk rises less steeply 
with the temperature, and the spectral shape is moderately distorted from that 
of the standard accretion disk. In this regime, radial temperature gradient of 
the disk has been found to be flatter than that of the standard accretion disk. 
The results of the apparently standard regime are suggestive of a slim disk (e.g., 
Abramowicz et al. 1988, Watarai et al. 2000) which is a solution predicted under 
high mass accretion rate. In the intermediate anomalous regime, the spectrum 
becomes much harder, and the disk inner radius derived using a simple disk 
model spectrum apparently varies significantly with time. These properties can 
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be explained as a result of significant thermal inverse Comptonization of the disk 
photons, as was found from GRO J1655 — 40 in its anomalous regime by Kubota, 
Makishima and Ebisawa (2001). 

Subject headings: accretion, accretion disks — black hole physics — stars: individual 
(XTE J 1550 - 564) —X-rays: stars 

1. Introduction 

In a close binary consisting of a mass accreting stellar-mass black hole and a mass do- 
nating normal star, the accreting matter releases its gravitational energy as X-ray radiation. 
When the mass accretion rate M is high, such a black hole binary is usually found in a 
so-called high/soft state, of which X-ray spectrum is characterized by a very soft component 
accompanied by a power-law tail. As described, e.g., by Makishima et al. (1986), the soft 
spectral component is interpreted as thermal emission from an optically thick accretion disk 
around the black hole, as it can be well reproduced by a multi-color disk model (MCD model; 
Mitsuda et al. 1984). This model approximates a spectrum from a standard accretion disk 
(Shakura & Sunyaev 1973). The MCD model has two spectral parameters; the maximum 
disk color temperature, T in , and the apparent disk inner radius, r in ; the latter can be related 
to the true inner radius, R[ n , via a correction for spectral hardening (Shimura & Takahara 
1997) and a boundary condition (Kubota et al. 1998). As the disk luminosity changes signif- 
icantly, the value of R m is usually observed to remain constant at the innermost Keplerian 
orbit for the black hole, QR g where R g = GM/c 2 is the gravitational radius (e.g., Ebisawa 
et al. 1993). 

Although this "standard picture" is successful for many of the high/soft-state black hole 
binaries (e.g., Tanaka & Lewin 1995; McClintock & Remillard 2003), it has been pointed out 
theoretically that the standard disk can be stable only over a limited range of M. Actually, 
a series of new solutions to the accretion flow have been discovered, including a slim disk 
solution, which takes advective cooling into account (Abramowicz et al. 1988; Watarai et 
al. 2000). In addition, a "very high state" has been observationally found as a derivative 
state from the soft state (Miyamoto et al. 1991; van der Klis 1994). Characterized by an 
enhanced hard component and significant variations in both T in and r in , the very high state 
is regarded as a possible violation of the simple-minded standard-disk picture. 

A clue to this problem has been recently given by Kubota, Makishima, & Ebisawa 
(2001; hereafter Paper I) from an analysis of the multiple RXTE/PCA data of the black 
hole transient with superluminal jet, GRO J1655 — 40. While the source behavior was 
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described adequately by the standard-disk picture over some period (called standard regime) 
of the entire PCA data span, the other period was characterized by the previously suggested 
deviation from such a standard behavior (called anomalous regime). The enhanced hard 
X-ray spectrum in the anomalous regime has been interpreted successfully as a result of 
significant inverse Compton scattering of the disk photons by some high energy electrons. 
The inner radius of the underlying optically thick disk has been found to be kept constant, 
when the effect of the Comptonization is taken into account. This result has been reinforced 
by Kobayashi et al. (2003). 

In order to reinforce the view obtained from GRO J1655 — 40, and to deepen our 
understanding of the physics of accretion under high values of M, the RXTE data of the X- 
ray transient XTE J1550 — 564 is analyzed in this paper. This transient source was discovered 
on 1998 September 7 by RXTE/ ASM and CGRO/B ATSE (Wilson et al. 1998, Smith 1998), 
and now confirmed as a superluminal jet source (Hannikainen et al. 2001). Figure 1 shows 
a lightcurve of this source, obtained with the RXTE/ ASM. As indicated with down-arrows 
in this figure, this outburst was continuously monitored by the RXTE pointing observations 
through 1999 May 20. Optical observations have established that the system consists of a 
late type sub giant (G8IV to K4III) and a black hole, the latter having a dynamical mass 
of Mx = 8.4-11.2 M Q (Orosz et al. 2002). The binary inclination angle and the distance to 
the source are estimated to be i = 67°-75° and D = 2.8-7.6 kpc, respectively. In this paper, 
i = 70° is used as a crude estimate, and its distance is denoted as D = D 5 ■ 5 kpc. 

In §2, the observation and data reduction are briefly described. In §3, the PCA spectra 
are analyzed using the canonical MCD plus power-law model, leading to the identification 
of three characteristic regimes; standard regime, anomalous regime, and apparently standard 
regime, in the increasing order of the luminosity. It is confirmed in §4 that the spectra in 
the anomalous regime can be well explained by the strong inverse Compton scattering, as 
found in Paper I. In §5, the apparently standard regime spectra are characterized in terms 
of the radial temperature gradient of the optically-thick disk, with a conclusion that a slim 
disk is probably realized. 

2. Observation and data reduction 

As indicated with down-arrows in Fig. 1, the entire outburst was covered by 184 point- 
ing observations with RXTE; the obtained data was analyzed in a standard way by Sobczak 
et al. (2000). By analyzing the outburst rise phase (~40 days in Fig .1) covered by 14 
pointing observations, Wilson & Done (2001) reported that this source experienced a sig- 
nificant spectral evolution. Figure 1 in their paper shows that the power-law photon index 
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describing the PCA spectra below ~ 20 keV switched from 1.4-1.7 to 2-3, which are typical 
of the low/hard state (e.g., Tanaka 1997), and the high/soft or very high state (e.g., Grove 
et al. 1998), respectively. Meanwhile, the 20-100 keV portion of the spectrum maintained a 
convex shape. Therefore, the source is inferred to have made a spectral transition from the 
low/hard state into the very high state (Wilson & Done 2001), or into the anomalous regime 
defined in Paper I, rather than into the canonical high/soft state wherein the spectrum above 
20 keV would not exhibit any cutoff. 

These published results indicate that the spectra acquired during the first ~ 40 days 
of the outburst is possibly related to development of the optically thick accretion disk; in 
response to a sudden increase in M, an optically thick accretion disk developed inward and 
reached the last stable orbit. Also the data after day 225 (1999 April 20) shows characteristics 
of the low/hard state. Therefore, the present paper focuses on the data taken from 128 PCA 
pointing observations between day 43 (1998 October 20) and 224 (1999 April 19), which are 
thought to represent relatively steady states with high M. For the systematic analyses of 
all the data set, the HEXTE data is not analyzed here because of poor statistics in some 
pointings. 

Following the standard procedure for bright sources, good PCA data was selected and 
processed. The data was excluded when the target elevation angle was less than 10° above the 
Earth's Limb, or when the actual pointing direction was more than 1'.2 away from the pointed 
direction. In particular, the data was discarded if it was acquired within 30 minutes after the 
spacecraft passage through South Atlantic Anomaly. The selected data from the individual 
proportional counter units was co-added and used for spectral analyses. The standard dead 
time correction procedure was applied to the data. The PCA background was estimated 
for each observation, using the software package pcabackest (version 2.1e), supplied by the 
RXTE Guest Observer's Facility at NASA/GSFC. Sometimes, the use of pcabackest resulted 
in an overestimate of the PCA background up to ~ 10%. Such a systematic over-estimate 
was corrected in the same way as in Paper I. That is, the on-source spectra were compared 
to the predicted model background spectra in the hardest energy band (>80 keV), where 
the signal flux is usually negligible. If necessary, the normalization factor of the background 
spectrum was changed. The PCA response matrix was made for each observation by utilizing 
the software package pcarsp (version 7.10). In order to take into account the calibration 
uncertainties, 1% systematic errors are added to the data. Over the 20-35 keV range, the 
systematic errors are increased to 10%, to cope with the response uncertainties associated 
with the Xe-K edge at ~ 30 keV. Although this could be an overestimate, the 20-35 keV 
data is utilized only in §4, and the results remain essentially unchanged even if it is reduced 
to 2 %. 
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3. Standard modeling of the PCA spectra 

3.1. Characterization of the observed spectra 

The 3-20 keV PCA spectra of XTE J1550 — 564 were analyzed by employing the canoni- 
cal MCD plus power-law model. The two constituent continuum components were subjected 
to common photoelectric absorption, with the column fixed at Nu = 1 x 10 22 cm -2 , which is 
reasonable for the source distance of 5 kpc and its Galactic position of (/, b) = (325°.9, — 1°.83). 
The following results are not affected by changing the value of Nu to 5 x 10 21 cm" 2 . The 
data requires additionally an absorption edge to the power-law component around 7-9 keV 
in terms of a smeared edge model (Ebisawa et al. 1994) and a narrow Gaussian line around 
6.5-6.7 keV for the Fe-K line. 

Figure 2 shows time histories of the best-fit model parameters, including the disk 
bolometric luminosity Ldisk (Mitsuda et al. 1984; Makishima et al. 1986), the 1-100 
keV power-law luminosity L pow calculated assuming an isotropic emission and their sum, 
L to t = £disk + £pow Along the time history, the entire observational span can be divided into 
the following eight periods (Periods 1-8 denoted in Fig. 2a) of relatively distinct properties, 
by mainly referring to the behavior of L pow and L disk . 

• Periods 1 and 7 are characterized by dominance of L pow and significant reductions in r in . 

• In Periods 2 and 8, L tot decreases keeping L pow rather low, while r in remains rather large 

and constant. 

• In Period 3, Ldisk gradually increases, but L pow remains negligible. The power-law com- 

ponent is too weak to constrain its photon index T. It is therefore fixed at T = 2.0. 

• In Periods 4 and 5, L pow is negligible (hence T being still fixed at 2.0), while L tot saturates 

at ~ 6 x 10 38 • L> 5 2 erg s -1 ; this is about ~ 40 • D 5 2 % of the Eddington limit, L E , for 
a black hole of 10 M . 

• Period 6 may be intermediate between Periods 5 and 7. 

Typical spectra representing these periods are shown in Fig. 3, and their best-fit parameters 
are given in Table 1. The fits are acceptable for the data in Periods 1, 2, and 6-8, but 
sometimes unacceptable in Periods 3-5 (Fig. 2e). This problem is considered again in §5. 
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3.2. Observed three spectral regimes 

To examine each Period for the validity of the standard picture, Fig. 4 shows several 
scatter plots between the spectral parameters. The data points of Periods 2, 3, and 8 
are thus confirmed to satisfy the standard picture, because r in is therein kept constant at 
~ 60 ■ D 5 km while T in changes significantly, ~ 0.5-1 keV (Fig. 2, Fig. 4a). By taking into 
account a correction factor for the boundary condition of £ = 0.412 (Kubota et al. 1998) 
and a color hardening factor of k — 1.7 (Shimura & Takahara 1997), the true inner radius 
is estimated as Ri n = r in • £ • k 2 ~ 71 • D 5 km. This value seems to be slightly smaller than 
6R g for a 10 M Q black hole, 90 km, though the source distance is not well constrained. 

The data points in Period 1 and 7 clearly violate the standard picture. Even though 
the fits are acceptable (Fig. 2e), r in is observed to change significantly over 30-50-_D 5 km 
(Fig. 2c), and T- ni shows strong positive deviation from their long-term trend (Fig. 2b). In 
these Periods, the hard component dominates the MCD component (Fig. 2a), and values of 
T are largest among these eight Periods (Fig. 2b). All these properties make Periods 1 and 
7 reminiscent of the anomalous regime of GRO J1655 — 40 (Paper I). 

The classification of Period 4-6 is somewhat ambiguous. The hard emission is negligible 
(Fig. 2a), and the spectral shape is similar to those in the standard regime (Fig. 3). However, 
as is clear from Fig. 2a, the data suffers from a strong saturation in Ldisk- In addition, a slight 
increase of T in under the saturation in Ldisk gives rise to a slight decrease in r in (Fig. 4a). 
In this sense, Periods 4-6 are called the apparently standard regime in this paper. Thus, 
XTE J1550 — 564 exhibits three characteristic spectral regimes, the standard (Periods 2, 3, 
8), apparently standard (4, 5, 6), and anomalous (1, 7) regimes. 

The Ldisk-^pow diagram, presented in Fig. 4c, can be used to understand the spectral 
evolution. The entire span of this source is found to show a clockwise loop; starting from a 
highest-L pow point (Period 1; anomalous regime), the source moves to the left along a low- 
L disk branch (Period 2; standard regime), then increases in L disk with L pow kept low (Period 
3; standard regime), and reaches a ceiling at L disk , ~ 6 x 10 38 • D 5 2 erg s _1 . It then moves to 
the right (Periods 4, 5 and 6; apparently standard regime). Finally, in Period 7 (anomalous 
regime), it returns to nearly the same position as Period 1, with Period 8 (standard regime) 
being a simple repetition of Period 2. Thus, the source fortunately exhibited a complete one 
cycle. 

Figure 5a shows a Ldisk-^ln diagram, which is useful to examine the validity of the 
standard picture against luminosity. A simple relation of Ldisk oc T m A means the constancy 
of r; n (see also Fig 4a), and hence the goodness of the standard picture. In this diagram, the 
data points in the standard regime indeed satisfy this relation, while those in the anomalous 
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regime deviate significantly. The data points in the apparently standard regime, clustered at 
the uppermost end of the diagram, deviate weakly from the standard Ldisk ^ln 4 relation, 
exhibiting a flatter dependence of Ldisk on T in as L disk oc T in 2 . 

It is useful to note here that the overall spectral shape of the MCD model well describes 
emission from a standard disk, even though it ignores the inner boundary condition of the 
disk. As long as the PCA spectra of T in = 0.5-2 keV are concerned, the values of T in and 
r in (with the corrections mentioned before) agree with 4-5 % with those obtained by a more 
accurate model, e.g., a cliskpn model (Gierlihski et al. 1999; Gierlihski & Done 2003) in 
xspec. Therefore, the characteristics shown in Fig. 4-5, are not due to incompleteness of the 
MCD model but mean that the anomalous regime and the apparently standard regime are 
intrinsically different from the standard regime. 

4. Reanalyses of the anomalous regime data 
— confirmation of the strong inverse Compton scattering — 

Now that Period 1 and 7 have been inferred to be the anomalous regime, these data can 
be reanalyzed employing the concept of the strong disk Comptonization which successfully 
explains the anomalous regime of GRO J1655 — 40 (Paper I). In the case of GRO J1655 — 40, 
the hard component in the anomalous regime has been found to be different in behavior from 
the power-law tail of the standard regime, in several points, including negative correlation 
between L pow and Ldisk, and systematically higher values of T. It is hence argued in Paper I 
that in the anomalous regime the spectrum is contributed significantly by a third spectral 
component, which is harder than the MCD emission but softer than the power-law tail in 
the standard regime. The strong anti-correlation between L pow and Ldisk has been taken 
for evidence that this third component strongly and negatively correlates with the MCD 
component. It is therefore natural to assume that a fraction of the photons from the optically 
thick accretion disk are converted into the third spectral component, most probably through 
inverse Compton scattering by high energy electrons which may reside somewhere around 
the disk. 

Following Paper I, the PCA spectra in the anomalous regime of XTE J1550 — 564 
are re-fitted with a three-component model, obtained by adding a Comptonized component 
to the original two component model. In this paper, a thermal Comptonization model 
(thcomp; Zycki, Done, & Smith 1999 ) is utilized to reproduce the Comptonized component, 
instead of the Comptonized blackbody model (compbb; Nishimura, Mitsuda, & Itoh 1986) 
which was used in Paper I. The thcomp model is based on a solution of the Kompaneets 
equation (Lightman & Zdziarski 1987). It well describes thermal cutoff around the electron 
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temperature of the plasma, T e , and allows the MCD model to be used as a seed photon 
spectrum. In contrast, the compbb model is useable only energies below kT e , and it assumes 
a single temperature black body emission as a seed photon spectrum. The thcomp model 
has actually improved the fit goodness compared to the compbb model because of difference 
in the seed photon spectra. The main results of the present paper are however independent 
of such modelings, and same results were obtained by using the compbb model. In this 
subsection, the 3-50 keV data is used to better constrain the wide-band spectral shape. 

For the spectral fitting, many of the model parameters were fixed to default values 
after Paper I. Namely, the maximum color temperature of the seed photons was tied to T[ n , 
to reproduce the situation whereby a part of the original MCD photons are up-scattered. 
Furthermore, T e was fixed at a representative value of 20 keV since it could not be constrained 
and is consistent with that obtained during the first 14 observations by Wilson & Done 
(2001). Neither a reflection component nor relativistic smearing was added to the model. 
Hence the free parameters of the thcomp model are two; thcomp photon index r tnc which 
expresses the spectral shape below kT e , and its normalization. Moreover, a value of T of the 
original power-law component was fixed at 2.0. As a result, the additional number of the 
free parameters for spectral fittings is reduced to only one. 

Figure 6 shows the same anomalous-regime spectrum as presented in Fig. 3a, but fitted 
with the three-component model over the expanded energy range. A solid line represents 
the additional thcomp component. The fit is acceptable, and the result implies that the 
dominant hard spectral component is mostly produced by the strong disk Comptonization. 
The re-estimated values of T in for all the data in the anomalous regime are plotted again in 
Fig. 2b with open triangles. By considering the disk Comptonization, the highly deviated 
data points in terms of T in have thus settled back to a smooth long-term trend as was 
already found in GRO J1655 — 40 (Paper I). The luminosity is also re-estimated in Fig. 5b 
as -^disk + L thc , where L thc is the estimated 0.01-100 keV thcomp luminosity, assuming an 
isotropic emission. Thus, Ldisk + -^thc plotted against the revised T; n approximately recovers 
the standard Ldisk oc T in 4 relation for optically-thick accretion disks. 

The value of r in is difficult to be estimated precisely under the strong Comptonization, 
because of uncertainties of both geometry and optical depth, r cs , of the cloud. In the present 
paper, it is approximately calculated by referring to observed photon flux and a re-estimated 
value of Ti n , on an assumption that few photons are scattered back into the optically thick 
disk and hence the number of the observed photons is conserved through inverse Compton 
scatterings. Details of this procedure are described in Appendix-A. The re-estimated values 
of ri n , plotted in Fig. 2c with open triangles, now appear to remain almost stable at ~ 
60 • D 5 km. Consequently, the optically thick accretion disk can also be considered to remain 
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relatively stable, even when a significant fraction of the MCD photons are Comptonized. 
Therefore, the picture of disk Comptonization is reconfirmed in the anomalous regime, as 
suggested in GRO J1655 - 40 (Paper I; Kobayashi et al. 2003). 

Table 2 summarizes the best-fit parameters associated with the exemplified spectra in 
the anomalous regime. The value of r t h c ~ 2.8 implies a ^/-parameter of ~ 0.2, or t cs ~ 1.8 
calculated via a following formula (e.g., Sunyaev & Titarchuk 1980) assuming T e = 20 keV; 



' O OK J 1 K 

■ + (T c /511keV) • {(r thc + 0.5) 2 - 2.25} '° 



The derived parameters are similar to those of GRO J1655 — 40. The smallness of y and r cs 
is consistent with the comparatively small fraction of L thc relative to Ldisk + L thc , typically 
< 0.5. It is also consistent with the assumption made in re-evaluating r in in Appendix-A 
and Fig. 5b, that the mean number of scattering is not too large and the fractional energy 
change is small. 

From these results, the anomalous behavior of XTE J1550 — 564 observed in Periods 
1 and 7 can be identified with that of GRO J1655 — 40 in the anomalous regime: the 
optically-thick standard accretion disk is present, but the Comptonization converts a signif- 
icant fraction of its emission into the hard component. 



5. Reanalyses of the apparently standard regime data 

5.1. Properties of the apparently standard regime 

As shown in §3, the apparently standard regime corresponds to the most luminous phase 
of the outburst, and the data in this regime occupies the upper-right region of the L^isk-Tm 
diagram (Fig. 5b). In this regime, Tj n gradually changed keeping Ldisk almost constant at 
6 x 10 38 • -D5 2 erg s™ 1 (~ 0.4 Le). As a result, the data points deviate from the standard 
Ldisk oc T in 4 relation, as L disk oc T in 2 . The obtained values of r in are not constant but 
become smaller than those in the standard regime, exhibiting a weak correlation with T in as 
r in oc Tin™ 1 . Figure 4a clearly shows this behavior. 

As seen in Fig. 3, the spectra in the apparently standard regime show a dominant soft 
component accompanied by a very weak hard tail. Although these properties are similar to 
that of the standard soft state, the apparently standard regime is something different from 
the standard regime because of the inconstancy of r in (or moderate saturation of Ldisk) , and 
the absence of the Fe-K line feature that is usually found in the standard regime. Moreover, 
the canonical spectral model often failed to give acceptable fits to the data (Fig. 2e). As are 
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clearly seen in the residuals of Fig. 3d-e, the discrepancy between the data and the best-fit 
canonical model appears as a low-energy excess and a spectral hump around 13 keV. This 
could be a result of fixing T of the hard tail component at 2.0; accordingly, the fits were 
repeated by leaving T free to vary. Then, the fits became acceptable, but they required 
so large values of T (see Fig. 2d, Table 1), that the spectrum below ~ 5 keV is mostly 
accounted for by the power-law component rather than the MCD component; the observed 
low-energy excess is filled up artificially by the steep power-law. Such a fit could be physically 
inappropriate. 

The above results on the apparently standard regime suggest a subtle difference in the 
accretion disk configuration from the standard disk, in such a way that the softest end of the 
observed spectrum is more enhanced than is described by the MCD model. Such a change, in 
turn, may arise if, e.g., the radiative efficiency of the inner disk region becomes reduced and 
the radial temperature gradient flattens. In order to quantify this idea, in §5.2 a p-free disk 
model is constructed as a mathematically generalized function of the MCD model. The data 
in the apparently standard regime is then examined in §5.3 by utilizing this model function. 



5.2. Formalism of the p-free disk model 

The concept of the standard accretion disk assumes that the energy released by accretion 
is half stored in the Keplerian kinetic energy, and half radiated away as local blackbody 
emission. As a result, the spectrum of the standard disk can be described as a geometrically- 
weighted sum of multi-temperature blackbody components of which the local temperature 
depends as T(r) oc r~ 3 / 4 on the distance, r, from the central black hole. Therefore, any 
departure of the physical condition assumed by the standard disk picture will make the 
radial temperature gradient deviate from the canonical value of —3/4. Such a deviation will 
in turn cause a slight deformation of the radially-integrated X-ray spectra that are observed. 
In order to quantify this idea, the MCD model is generalized as seen below, after initial 
attempts by Mineshige et al. (1994, for GS 1124-68) and Hirano et al. (1995, for Cyg X-2). 

The main assumptions of the model developed here are that a disk local temperature 
is described by T(r) = T in ■ (r/r in ) _p , and that the disk locally emits a blackbody spectrum. 
Here, p is a dimension-less positive parameter introduced to generalize the MCD formalism, 
with p = 3/4 implying the MCD model. The spectrum from this model function can be 
written as 

<M-*?P 

with B(E,T) being a blackbody spectrum of temperature T. For convenience, hereafter 
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this mathematical model function is called p-free disk model. As p decreases, the spectrum 
becomes softer than the MCD spectrum of the same T in , because the radiation from outer 
parts of the accretion disk is emphasized as seen in equation (2). 

5.3. Spectral fitting with the p-free disk model 

The PCA spectra of Periods 3-5 (one segment of standard regime and two of apparently 
standard regime) have been re-fitted by replacing the MCD model component with the p- 
free disk model. Although an accurate determination of p is difficult when the power-law 
component is strong, the PCA spectra in Period 3-5 are fortunately free from this obstacle. 
The condition of spectral fitting is otherwise the same as in §3. That is, the absorption 
column iVn and T of the power-law component are fixed at 1 x 10 22 cm -2 and 2.0, respectively. 
The Gaussian line was not included. 

The time histories of the obtained p-free disk model parameters are given in Fig. 7. In 
Period 4 and 5 (apparently standard regime), the fit goodness has been significantly improved 
by allowing p free. Figure 8 shows the same apparently standard regime spectrum originally 
presented in Fig. 2e, fitted this time with the p-free disk model, and Table 3 shows the 
examples of the p-free model parameters. 

In Fig. 9, the best fit values of p are plotted against T- m ; here, instead of the values of T m 
obtained by the p-free disk model, those by the MCD model are employed, in order to avoid 
any systematic coupling between p and T in . The dependence of p on T in thus changes at ~1 
keV. As a function of T in , p increases up to ~ 1 keV (Period 3, the standard regime), beyond 
which it decreases abruptly and the correlation turns negative (Period 4-5, the apparently 
standard regime). To examine these results for various systematic effects, the p-free disk fits 
have been extensively repeated by changing the fitting conditions; fixing JVh to 5 x 10 21 cm -2 
instead of 1 x 10 22 cm -2 , fixing T to 2.2 instead of 2.0, including a Gaussian line of which 
the central energy is constrained to 6.2-6.9 keV. These different conditions slightly affected 
the absolute values of p, but did not affect the characteristic p vs. T in behavior of Fig. 9a. 

Thus, the value of p has been found to deviate from 3/4 as the source enters deep into 
the apparently standard regime. However, a still larger excursion of p is observed during the 
standard regime, where p = 3/4 should be obtained. Therefore, the changes of p in Fig. 9 
could partially or entirely be due to some artifacts, and p could deviate from 3/4, even if the 
standard accretion disk is realized. This could actually happens, because the MCD model is 
a mere approximation of the exact standard-disk solution; the actual temperature gradient 
of a standard accretion disk must be flatter than —3/4 near the innermost disk edge, where 
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the temperature will approach zero. As T in decreases, the limited PCA band pass will sample 
preferentially the emission from inner disk regions, thus making p different from 3/4. 

In order to address the above issue, a theoretical approach was first attempted. That 
is, a number of 3-20 keV PCA spectra were simulated, using two theoretical model spectra 
in the xspec which are known to be more accurate than the MCD model. One is the diskpn 
model, which is based on the Shakura-Sunyaev solution in a pseudo-Newtonian potential; it 
properly takes the inner boundary condition into account. The other is so-called a GRAD 
model (Hanawa 1989; Ebisawa, Mitsuda, & Hanawa 1991), which considers full general 
relativistic effects for a Schwarzshild metric. Then, the simulated spectra were fitted with 
the p-free disk model. As a result, both diskpn and GRAD models gave p = 0.72-0.75 
as long as T in is in the range 1.2-1.5 keV. Furthermore, p was found to change (over 0.5- 
1.0) for lower input disk temperature of T in = 0.8-1.5 keV. However, the p-T- m relation 
turned out to be quite different between the two input models, and to depend significantly 
on the inclination angle when the GRAD model is used. In addition to these theoretical 
complications, neither the diskpn model nor the GRAD model can be considered accurate 
yet, since they (as well as the MCD model) neglect possible spectral deviation from a pure 
local blackbody. Accordingly, this approach has been concluded unrealistic. 

In this paper, instead, the PCA-determined values of p for a standard disk has been 
calibrated empirically, by using actual spectral data of a prototypical black hole binary, 
LMC X-3, in the standard regime. This black hole binary perfectly satisfies the standard 
picture up to ~ Le for a black hole of 5-7.2 M Q (Paper I, Wilms et al. 2001), and its system 
inclination angle, 65°-69°, is similar to that of XTE J1550 — 564. The p-free disk model 
was applied to the PCA spectra of LMC X-3 obtained by 128 pointed RXTE observations 
from 1996 February to 1999 January (Kubota 2001). The best fit values of p from LMC X- 
3 are plotted against T in in Fig. 9b. As expected, a positive correlation between p and 
T in artificially appears even for this prototypical "standard- disk" object. Therefore, the 
behavior of XTE J 1550 — 564 (in Fig. 9a) for T m < IkeV (Period 3, the standard regime) can 
be understood as an artifact. In contrast, the behavior of XTE J1550 — 564 for T m > 1 keV 
is clearly distinct from that of LMC X-3. The conclusion is that the temperature gradient 
of XTE J1550 — 564 in Period 4 and 5 becomes intrinsically smaller than in the standard 
regime, and hence the accretion disk in the apparently standard regime in reality deviates 
from the standard picture. 
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6. Discussion 

6.1. Overall picture from the observation 

Through the detailed analysis of the PCA data of XTE J1550 — 564, the three spectral 
regimes have been identified. Their relation on the L disk -T in plane is illustrated schematically 
in Fig. 10, while their properties can be summarized as follows. 

1. When Ldisk is well below a certain critical upper-limit luminosity, L c ~ 6 x 10 38 • 
D 5 2 erg s _1 (~ 0.4 • D 5 2 Le), the spectral behavior can be explained by the standard- 
disk picture. This is the standard regime. 

2. When Ldisk hits L c , it is moderately saturated as Ldisk oc T in 2 , and hence r in shows 
a weak correlation to Xj n T in 1 - Although the spectrum, consisting of a 
dominant soft component and a weak hard tail, resembles that in the standard regime, 
the radial temperature gradient in the disk (represented by p) becomes flatter than in 
the standard disk. This is the apparently standard regime (§5). 

3. At an intermediate case (Ldisk ~ L c ), the spectral hard component dominates, and 
the apparent inner-disk radius is no longer constant. This is the anomalous regime 
(§4). These effects can be explained by a sudden increase of the disk Comptonization, 
while the underlying disk itself remains in the standard state and r in is kept constant. 

4. The source evolves from the standard to apparently standard regimes, then to the 
anomalous regime, and returns again to the standard regime. 

The anomalous regime is naturally identified with that found in GRO J1655 — 40, and the 
scenario of Comptonization suggested in Paper I successfully applies to XTE J 1550 — 564 
as well. It has been confirmed that the violent variation in r- m in the anomalous regime is 
apparently caused by strong disk Comptonization, with the underlying optically thick disk 
extending down to the last stable orbit like in the standard regime. The apparently standard 
regime is possibly the same as Period 1 of GRO J1655 — 40. Because intensity variation 
in Period 1 of the source was small, Paper I did not discuss the source behavior in that 
period. However, GRO J1655 — 40 resides in this period in the upper right corner on the 
£disk-Ti n plane (see Paper I), and its spectra consist of a dominant disk component and a 
very weak hard tail. These properties are basically the same as those of XTE J1550 — 564 
in the apparently standard regime. 
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6.2. Comparison with theoretical predictions 

As is well known, theoretical solutions to the steady-state accretion flow form an S- 
shaped locus on the plane of M vs. the surface density of the disk (e.g., Abramowicz et 
al. 1988, Chen & Taam 1993, Kato, Fukue, & Mineshige 1998). The locus involves to 
thermally stable branches; the standard Shakura-Sunyaev accretion disk solution, and the 
slim-disk solution, realized when M is relatively low and very high, respectively. Evidently, 
the standard regime can be identified with the standard Shakura-Sunyaev solution. 

The slim disk solution takes into account the effect of the advection, in addition to the 
viscous heating and the radiative cooling. This effect becomes important when M is high 
and hence the luminosity is close to Le- Under this condition, any increase in M would be 
balanced by an increase in the advective transport, accompanied by little increase in L disk ; 
this agrees with the observed mild saturation in Ldisk observed in the apparently standard 
regime. Watarai et al. (2000) simulated many slim disk spectra, fitted them with the MCD 
model, and derived an empirical relation between r in and the apparent T in as r in oc Tj n _1 (or 
-^disk Tm 2 ). This is exactly what has been observed in Fig. 5a. Furthermore, Watarai et 
al. (2000) showed that the temperature gradient becomes flatter as the advective cooling 
becomes important because of a progressive suppression of disk emissivity. In the extreme 
case, p can reduce to 0.5. Thus, the overall source behavior in the apparently standard regime 
agrees very well with the prediction by the slim disk model. Of course, this state assignment 
is still tentative, because the slim disk solution still neglects important effects such as general 
relativity, magnetic field, and photon trapping (Ohsuga et al. 2002). Other solutions would 
have to be considered as well. 

In addition to the two stable branchs, a thermally and secularly unstable branch is 
known to exist between them. This branch is recognized as a negative slope on the S- 
shape sequence. By comparing the obtained Ldisk-^ln diagram to the S-shape sequence, the 
anomalous regime may have some relation to the unstable branch of the sequence. In other 
words, the instability of the standard disk may cause the anomalous regime. Interestingly, 
the quasi-periodic oscillations (QPOs) are observed preferentially in the anomalous regime, 
in GRO J1655 - 40 (Remillard et al. 1999) and XTE J1550 - 564 (Remillard et al. 2002). 
Therefore, the QPO is likely to relate to the existence of the Compton cloud. 

The observed three distinctive regimes are likely to reflect the change of the accretion 
disk structure from the standard accretion disk to other solutions, as the radiative cooling 
becomes progressively inefficient. These results hence provide, at least potentially, one of 
the first observational accounts of the long predicted S'-shaped sequence. 
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A. A formula to re-estimate r in under Comptonization 

Under the presence of strong disk Comptonization, the apparent disk inner radius r in 
may be calculated as 

FL k + F L • 2COS. = 0.0165 . (jgj^) ■ (j^)' P^ons r' cm"» , (Al) 

where F^ isk and F^ hc are 0.01-100 keV photon flux from the direct disk component and the 
Comptonized component, respectively. Here, T in refers to the disk temperature obtained by 
considering the inverse Compton scattering as is done in § 4. The first parentheses is just the 
same as the normalization factor of the MCD model in the xspec. This formula is based on 
an assumption that there are few photons which are injected again into the optically thick 
cool disk (i.e., optical depth of the cloud is not very large). It can be derived through the 
following steps. 

1. Via flux command of the xspec, 0.0165 photons s _1 cm~ 2 in the range of 0.01-100 keV 
is obtained for the MCD spectrum of which T in and normalization are 1 keV and 1, 
respectively. 

2. According to the Stefan-Boltzmann's law, photon flux from the original optically thick 
disk emission is in proportion to rf n • . 

3. Isotropic emission is assumed for Comptonized photons, and thus F^ hc is multiplied by 
2 cost. 
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Table 1: The best-fit parameters of XTE J1550 - 564 with the MCD fit 



Per. 


date a 


Tin [keV] 


r in [km] b 


r 


Tdisk ° 


T c 

^pow 


E c [keV] 


eqw [eV] 




1 


43 


1.05 ± 0.04 


36^ 


2.47 ±0.05 


2.07 


2.18 


6 6 +u ' 3 


45 ±22 


15.3/37 


2 


71 


0.67 ±0.01 


60 ±3 


1 Q7+0.07 
— 0.06 


0.92 


0.062 


6 5+ 02 


84 ±28 


21.2/37 


3 


91 


0.682 ± 0.008 


60 ±2 


(2.0 fixed) 


1.01 


0.027 


fi O+0.3 
D.J_ 2 


80 ±30 


24.9/38 






0.6861°;^ 




n 9 + O.I 
z - -0.2 


1.02 


0.028 


6.2 ±0.3 


68t1o" 


23.1/37 


4 


125 


1.117 ±0.004 


53.7 ±0.6 


(2.0 fixed) 


5.80 


0.022 




< 10.3 d 


94.0/38 






1 105+ ' 006 
±.±uo_ 004 


54 ±1 


o 7 +0.5 
J -'-0.3 


5.73 


0.369 


6 l +0 ' 5 

D - ± -0.4 


22 ± 17 


18.1/37 


5 


143 


1.141 ± 0.004 


52.1 ± 0.6 


(2.0 fixed) 


5.93 


0.031 




< 9.2 d 


92.8/38 






-, ioQ+0.007 
1.100_q 00 g 


51 ± 1 


3.9 ±0.3 


5.61 


0.834 


6 1+ ' 5 

"^-o.e 


18+?° 


21.4/37 


6 


172 


1.06 ±0.02 


57 ±2 


9 91+O.II 
z - zi -0.08 


5.27 


0.75 


6.7 ± 0.3 


46±?I 


9.91/37 


7 


182 


1.10 ± 0.05 


31 ±2 


9 cr+0.08 


1.90 


2.85 


6 6+ - 4 


35 ± 19 


7.5/37 


8 


192 


0.78 ± 0.02 


60 ±3 


2 oi+° 11 


1.68 


0.19 


6.6 ±0.2 


63 ±24 


17.9/37 



NOTE. — Errors represent 90% confidence limits. The PCA 3-20 keV spectra are used. For all the spectral fitting, 
the smedge and a Gaussian with fixed a = 0.1 keV are included, and ./Vjj is fixed at 1 X 10 22 cm -2 . 



a Days since 1998 September 7. 

b Apparent inner radii under assumptions of i = 70° and D = 5 kpc. 

c Tdisk and Lpow represent the bolometric disk luminosity, 1-100 keV power-law luminosity, 

respectively, in the unit of 10 38 ■ -D5 2 erg s -1 . As for Ldisk> * = 70° is assumed. 
d Upper limit of Gaussian line with fixed E c = 6.4 keV. 



Table 2: The best-fit parameters with the thcomp fit 



Per. 


date 


T in [keV] 


r in [km] a 


r t hc 


-^disk b 


T b 


Lthc b 


X> 


1 


43 


96 +U ' 04 


53+ b 3 


2.8+11 


2.29 


0.70 


1.27 


33.6/74 


7 


182 


n 9fi+ 05 

u - yD -0.06 


55±| 


2.8 ±0.1 


2.37 


0.69 


1.73 


34.5/74 



NOTE.— The PCA 3-50 keV spectra are used. 

a Apparent inner radii re-estimated by equation (2). 

b ^diski L pow , and L t h c represent the bolometric disk luminosity, 1-100 keV power- law 
luminosity and bolometric (0.01-100 keV) thcomp luminosity, respectively, all in the 
unit of 10 38 ■ D 5 2 erg s 1 . Isotropic emission is assumed for L t hc an d ipow, and 



i = 70° is assumed for Ldisk- 
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Table 3: The best-fit parameters with the p-free disk fit 



Per. date T in [keV] p E c cqw [kcV] \ 2 /v 



3 


91 


76+ ,UHi 

u - '"-0.02 


41+ u - 04 

u '^ 1 -0.05 






25.4/39 






0.73 ±0.05 




u -°-o.i 


6 3+ ' 3 


48 


22.0/37 


4 


125 


1 16+ - 01 


59+ 04 






54.8/39 






-i 17 +0.()2 
1 - 1 '-0.01 


0.56 ± 0.03 


5 9+ - 2 


45 


38.9/37 


5 


143 


1 1S+0.O2 
1-1»_0.01 


0.60 ± 0.03 






57.6/39 






1.20 ±0.02 


0.57 ±0.03 


6 0+ ' 2 


31 


47.7/37 



NOTE. — The PCA 3-20 keV spectra are used. The power-law photon index is fixed at 2.0. 
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Days since 1998 September 7 

Fig. 1.— The 1.5-12 keV lightcurve of XTE J1550 - 564, obtained with the RXTE/ ASM. 
The pointing observations are indicated with down-arrows. The horizontal arrow represents 
the period studied in the present paper. 
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Fig. 2. — Evolution of the spectral parameters of XTE J1550— 564. (a) Time histories of L disk 
(filled circles), L pow (open circles), and L tot (a solid line), all in the unit of 10 38 • D 5 2 erg s _1 . 
As for the disk luminosity, % = 70° is assumed. The eight characteristic periods are indicated 
at the top. (b)-(e) Those of T in , r in , T, and x 2 /dof, respectively. Open triangles in panel 
(b), (c), and (e) are obtained by incorporating the disk-Comptonization. Crosses in panel 
(d) and (e) for Period 3-5 show the results from the free-f fitting. 
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Fig. 3. — Typical PCA spectra of XTE J1550 — 564 in Periods 1-8, compared with the 
predictions of the best-fit MCD plus power-law model. Background has been subtracted, 
but the instrumental response is not removed. The lower panel of each spectrum shows fit 
residuals. 
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Fig. 4. — Several scatter plots among the spectral parameters. Panels (a), (b) and (c) shows 
r in against T in , T against L pow , and Ldisk against L pow , respectively. In panels (a) and (c), 
Period 1-8 are specified by eight kinds of symbols presented in panel (a); data points of 
standard, anomalous and apparently standard regimes are shown with filled symbols, crosses, 
and grey open ones, respectively. In panel (b), the data for Periods 3-5 are excluded because 
of the weakness of the power-law tail. 
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Tin[keV] Tin [keV] 

Fig. 5. — (a) The calculated Ldisk plotted against the observed T in . Periods 1-8 are specified 
by the same eight kinds of symbols as Fig. 4. The solid and dotted lines represent the 
Ldisk oc T in 4 and Ldisk oc T in 2 relations, respectively. The source distance and inclination are 
assumed to be 5 kpc and 60°, respectively, (b) Same as panel (a), but the data points in 
Period 1 and 7 (anomalous regime) are re-calculated considering the Comptonized component 
as Ldisk + -^thc- 
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Fig. 6. — The PCA spectrum in the anomalous regime (the same as Fig. 3a), fitted with the 
three-component model incorporating the Comptonized component (the medium hardness 
one shown with a solid line). 
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Days since 1998 September 7 

Fig. 7. — Time histories of the best-fit parameters of the p-free disk model (diamond). For 
comparison, the MCD parameters are also plotted with filled circles in the top and bottom 
panels. In the middle panel, p = 0.75 is indicated as a dotted line. 
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Fig. 8. — The PCA spectrum in the apparently standard regime (the same as Fig. 2e), fitted 
with the p-free disk model plus power-law (r = 2.0 fixed). 
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(a)XTE J1550-564 (b)LMC X-3 
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Tin [keV] Tin [keV] 

Fig. 9. — The best fit values of p, plotted against T m determined with the MCD model. 
The PCA results for XTE J1550 — 564 in Period 3-5 (panel a) are compared with those of 
TMC X-3 (panel b). Symbols for XTE J1550 - 564 are the same as in Fig. 3. In panel (b), 
data points with large error bars (Ap > 0.5) are not shown. As obtained in the MCD fit to 
the data of LMC X-3 (Kubota 2001), the values of JV H and T are fixed at 7.5 x 10 20 cm" 2 
and T = 2.5, respectively. 




Fig. 10. — A schematic classification of the three spectral regimes on the L^^-T^ diagram. 
Thick solid and dashed lines show the source behavior obtained under the MCD plus power- 
law fit. 



